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Non-zero curvature in a waveguide leads to the appearance of an attractive quantum potential which 
crucially affects the dynamics in matter-wave circuits. Using methods of supersymmetric quantum 
mechanics, pairs of bent waveguides are found whose geometry-induced potentials share the same scattering 
properties. As a result, reflectionless waveguides, dual to the straight waveguide, are identified. Strictly 
isospectral waveguides are also found by modulating the depth of the trapping potential. Numerical 
simulations are used to demonstrate the efficiency of these approaches in tailoring and controlling 
curvature-induced quantum-mechanical effects. 

Waveguides with non-zero curvature are basic constituents of matter-wave circuits in atom chip techno- 
logy^'^, as well as its ion^ molecular^, and electron^ counterparts. Their relevance is further enhanced by 
the development of flexible techniques to create optical waveguides for ultracold gases. In this context, 
waveguide trapping potentials can be engineered by a variety of methods including the time-averaging painted 
potential technique^, the use of an intensity mask^'^, and holographic methods, in particular, digital holography^. 
Circular ring traps have attracted a considerable amount of attention^'^°"^^ most recently to study Josephson junction 
dynamics^^'^^. Other curved waveguides have also been engineered, such as a stadium-shaped potential trap^^'^^. 

The propagation of matter- waves in bent waveguides generally differs from that in straight waveguides due to the 
appearance of a purely attractive local quantum potential of geometrical origin^^"^^. Under tight-transverse confm- 
ment, the magnitude of this curvature-induced potential (CIP) is proportional to the square of the curvature of the 
waveguide, and affects both the single-particle and many-body physics of the confined matter-waves^^"^^. As a result, 
the scattering properties of a curved tight waveguide are modified, e.g., by the appearance of bound states^^'^^. 
Advances in the design of bent waveguides, in which curvature-induced effects are tailored and suppressed, are 
required for the miniaturization of matter-wave circuits. It is to this problem that we turn our attention. 

Results 

In this manuscript we design bent waveguides for matter-wave circuits free from spurious quantum mechanical 
effects associated with CIPs. Three novel ideas are presented: (i) Exploiting the interplay of geometry and 
supersymmetry in quantum mechanics, we relate pairs of waveguides whose CIPs are isospectral and share 
the same scattering properties, (ii) We then identify waveguides which are reflectionless for coherent matter- 
waves at all energies, (iii) Furthermore, we show that by tailoring the depth of the waveguide trap, it is possible to 
cancel the CIP, rendering the dynamics of the guided matter-waves equivalent to that in straight waveguides. 

Let us consider the dynamics of matter-waves confined in a tight waveguide whose axis follows the curve y, 
parametrized as a function of the arc length qi by the vector r = r{qi), with tangent t(qi). We start by recalling the 
fundamental theorem of curves which asserts that a curve is completely determined, up to its position in space, by 
its curvature k and torsion t^°. Indeed, the expressions k = K{qi) and t = T(qi) constitute the natural intrinsic 
equations of a curve. A parametrization of the curve can be obtained by integration of the Frenet-Serret equations 
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where n = ((it/ dqi) / k (provided di/ dqi ^ 0) and b = t x n are the 
principal normal and binormal unit vectors, and the curvature and 
torsion at the point of arc length qi are defined as K(qi) = 
\dt{qi)/dqi\ and T{qi) = —h-dh/dqi, respectively. Let {q2, qs) be 
the transverse local coordinates and consider a transverse confining 
potential (q2, qs) such that in the limit of tight confinement 
the particle is bounded to y. Under dimensional reduction, the 
purely- attractive CIP emerges^ ^"^^'^^ 



(2) 



This result is independent of and holds in particular under an 
isotropic transverse harmonic confinement Ua,^ = mw\ {ql + qf) /2 
with ground state width Gq = [h/(m<jj j_)Y^^^^'^^. The conditions for 
the dimensional reduction to be valid explicitly read 



(3) 



where primes denote derivatives with respect to qi. 

We next pose the problem of identifying pairs of waveguides with 
the same scattering properties, and engineering a waveguide which 
minimizes the effect of the CIP. Generally, direct integration of the 
Frenet-Serret equations is not possible. However, the quantum 
mechanical behavior of matter waves bounded to isometric curves 
with different torsion but the same curvature remains the same^^'^^ 
because ^ is independent of t. In addition, matter-wave circuits in 
atom chips and optical realizations of waveguides are often assoc- 
iated with curves y on a plane, for which t = 0. As a result, we focus 
on planar curves, given by the parametrization r{qi) = {x{qi),y{qi)). 
We shall return to the case of t 0 whenever the curve y exhibits 
multiple points in which the waveguide self-intersects. The CIP 
depends only on and remains invariant under the mapping k 
sgn{g(qi))K with an arbitrary real function g{x), a symmetry which 
we shall exploit to engineer the guiding potential. Provided t = 0, 
and K(qi) 0 for all qi it is always possible to integrate the Frenet- 
Serret equations, and to find the natural representation of the curve 
in terms of the arc length 



(4) 



In what follows we shall take (xo,yo) ^ without loss of generality, 
and extend Eq. (4) to < 0 to sample the full range of 

Supersymmetric partner waveguides. In the Witten model of 
supersymmetric quantum mechanics (SUSY QM)^^'^^, a pair of 
SUSY partner Hamiltonians is considered with the factorization. 



x{qi)=Xo-\- [ cos( K{s)ds\ds, 



where the annihilation and creation operators are defmed by 



+ ^{qi) and , 



+ 0(<2i), in terms of the 



v2m v2m 
superpotential ^(qi). The SUSY partner Hamiltonians can be explicitly 



written as H+ = 

h 



where V+ {qi) : = {qi ) + 
^'{qi) are partner potentials. We consider the case in which 



2 , , 8V2m^,, 



8m 



(6) 



(7) 



where we have used the relation between the derivative of the 
superpotential and the commutator of the creation and annihilation 
operators in the second line^^. The relation between k±, can be further 
developed by making reference to the ground state i/^o of H- satisfying 
0, in terms of which the superpotential reads 



2m lAo 



- . Using this expression in (6), it follows that 



K\{qi) = K^_{qi)+8 



^0 V 'Ao 



(8) 



According to the fundamental theorem of curves, the shape and length 
of a planar curve is completely determined by its (single-valued and 
continuous) curvature. It follows that the SUSY partner potentials V+ 
are associated with the family of curves {y+} (with curvatures k± 



satisfying kI^ = 



2m 



-0)2 + 



2m 



O'), that we shall refer to as SUSY 



n n 

partner curves. This set can be extended to include curves associated 
with a family of shape-invariant potentials'^, as discussed in Methods, 
or using higher-order SUSY QM'^. 

Why are waveguides along these curves of interest? The main 
physical feature of SUSY partner curves is that they exhibit the same 
scattering properties, a distinguishing feature directly inherited from 



Let y± be open waveguides with finite curvature as qi 



±00, so that ('^i-^ ± (X))^(^{qi-^ ± oo) =:^\, and consider 
the scattering states of momentum k and energy £ = h^k^/{2m), with 
reflection and transmission amplitudes R±{k) and T±{k), respect- 

ively. It follows that R_ (k) = ^ — .^/^ (^) and T_ (k) = 



-ihk'/\/2m 



ll/2 



ihk/ \/2m 
T+{k) where k= [2m{E-^^_)y 

K that is, the reflection as well as the transmis- 



h and k' ^ 



0_ — ihk/ \/2m 

sion probabilities are the same for SUSY partner curves. Further, the 
Hamiltonians H± associated with SUSY partner curves are isospec- 
tral, except for the lowest energy level of H_ with zero -energy, which 
is absent in the spectrum of H+. 

Design of reflectionless curves. CIPs are of attractive character and 
as a result can lead to quantum reflection'^"'^. The dynamics of a 
guided matter-wave on a bent waveguide is generally affected by the 
curvature. We next illustrate the power of the SUSY partner 
waveguides in designing reflectionless curves. An obvious instance 
where the CIP vanishes is that of an infinite straight waveguide, with 
K- = Q and superpotential = A tanh aqi with A > 0. This 
configuration is of relevance to guided atom lasers^°'^\ and we wish 
to mimic it in bent waveguides. SUSY QM allows us to find SUSY 
partners which are reflectionless. In this case, V+{qi) is given by the 



modified Poschl-Teller potential Vcurv('2i) = 
that the curvature of the SUSY y+ curve reads 



r v(v + i) 



2m cosh {aqi 



k:+ {qi) =2a^ v(v + l)secha<^i, 



, so 



(9) 



where v is a positive integer. Provided that the dimensional reduction 
is valid, the transmission probability for a waveguide with curvature 



flie SUSY partner potentials are both induced by curvature, i.e. ^^^^^^^^ ^' l^+l ^ 



h 

y±=^j^ = k\. It foUows that SUSY partner Hamiltonians 

8m - 

are associated with curves whose curvatures k± are related by 



with 



sinh(7i;/c/a) ^ ^ / — ^ ^ ^ 

Li= and k=\/2mE/n. Such a waveguide becomes 

sinTTV 

reflectionless for integer values of v. Different reflectionless 
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(a) (b) 




(c) (d) 




Figure 1 | Supersymmetric reflectionless waveguides, (a) Waveguide 
whose CIP is the reflectionless Poschl-Teller potential with v = 1. Its SUSY 
partner curve is the straight line y_ with zero curvature, (b) The multiple 
point can be removed to obtain a simple waveguide using the curvature 
mapping K{qi) ^ sgn.{qi)K{qi) under which the CIP remains invariant, 
(c) For higher curvature values, as in the waveguide associated with the n = 
2 Sukumar reflectionless potential shown here, integration of the Frenet- 
Serret equations leads to curves with several multiple points, (d) Such CIPs 
can be engineered in a non-planar waveguide, with non-zero torsion t. 

waveguides are plotted in Fig. 1, where it is shown that the number of 
multiple points increases with the magnitude of the curvature. A 
simple waveguide without junctions can then be engineered by 
exploiting the invariance of the CIP with respect to changes in the 
sign of the curvature, or by considering a nonzero torsion t 0, 
whose realization might be achieved using an extended version of the 
painted potential technique^. We emphasize that there are infinitely 
many instances of reflectionless waveguides. Let us illustrate this by 
adapting the algorithm developed by Sukumar^^ to CIPs. One can 
construct reflectionless waveguides supporting n bound states. Let us 

2 2 

frx the bound state energies to be £„= rj with rj > 

2m 

^l-i> • • - ^i- The symmetric reflectioneless curves are described 
by the equation 

Kliqi) =8dl\ndet Dn, (10) 



where [Dn\^=^n] ^^Q^^[r]fl^ +{-\y^^Q^^(^-nfl^'^. For a 

single bound state, one obtains = %f]\stch.f]iqi, closely related to 
the SUSY curves associated with (9). Figure 1 (lower panels) shows 
the reflectionless curves corresponding to a Sukumar potential 
supporting two bound states with r\i = \ and r\2 — 3/2, where 
multiple points in 1(c) are avoided by a non-zero torsion 
{T = 20q~^) in 1(d). The axis of the associated waveguide follows 
the curve (x(s), y{s), xs) with (squared) curvature k^{s) = 
— Sm^{s)/h^y torsion t and arc length qi = VI-\-t^s. At variance 
with (9), the relative angle between the asymptotes can be tuned by 
adjusting the value of f]2 relative to f/i, which will allow for the 
engineering of reflectionless bends through a range of desired 
angles. Further examples of reflectionless waveguides can be found 
by using the infinite family of reflectionless potentials discussed by 
Shabat^^ and Spiridonov^^. The reflectionless character of the SUSY 
waveguides becomes apparent in the dynamics of guided matter 
waves. Figure 2 shows an elongated Gaussian beam being guided 
in a bent waveguide with curvature given by (9). y+ is 
asymptotically flat for qi ^ ±oo. For a general non-integer value 
of V, the traveling beam is substantially reflected off the bent 
region. For integer v there exists a delocalized critical bound- state 
with zero energy and the waveguide becomes reflectionless for all 
scattering energies. However, the degree of bending increases with v. 
As a result, reflectionless waveguides provide a remarkable 
counterexample to the common expectation that the reflection 
probability increases with the degree of bending of the waveguide. 
In addition, the numerical simulations correspond to the 
propagation in a waveguide with finite transverse width, for which 
the explicit form of the curvature induced potentiaP^'^^'^^ is more 
complex than that in Eq. (2) used to design the reflectionless SUSY 
waveguide, and where excitations of the transverse waveguide modes 
are possible. The fact that despite the finite transverse width the 
waveguide remains reflectionless illustrates the robustness of its 
design against imperfections. We also note that the reflectivity of 
the Poschl-Teller potential changes only gradually as v departs 
from an integer value. 

Canceling out the geometry-induced potential. A variety of 
experimental techniques to design matter-wave circuits, such as 
the painted potential technique, offers an alternative way to 
control the design of y: the modulation of the potential depth of 
the waveguide. Consider two arbitrary isometric waveguides y and 



Figure 2 | Scattering dynamics in bent waveguides. Sequence of snapshots of the time- evolution of the density profile of a wavepacket along a planar bent 
waveguide with the curvature (9) and v = 1/2 (left), and v = 1 (right), as in Fig. 1(b). Generally, the wavepacket is split by the CIP into a transmitted 
and a reflected component. Despite the high degree of bending shown in the inset, whenever v is an integer, the waveguide becomes reflectionless 
and exhibits unit transmission probability for all energies of the impinging matter-wave beam. The color coding varies from white to red as the 
probability density increases. The dimensions of each waveguide image are 908cro X 47 Gq and the time interval between successive images is 1920/co j_ . The 
initial wavepacket has FWHM = 235(7o and momentum (l/32)mco^(7o, and a = 1/8. 
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s = 0.90 



s = 0.95 



e = 0.97 




Figure 3 | Canceling out the curvature-induced potential. Elliptical 
waveguide potentials of increasing eccentricity (top) and corresponding 
ground state densities (middle). Bottom row shows ground state densities 
when the CIP is compensated by modulating the depth of the trap. The 
dimensionless density profile n{qi)aQ is scaled up by a factor 10^, the 
perimeter of the ellipse is L = ISOcto and the plotted area is SOctq X SOctq. 

y (both either open or closed, and without multiple points), with CIPs 
^y{q\) and ^y{qi), respectively. Under the consistency conditions 
(3), it is then possible to make y isospectral to y by modulating the 



depth of the waveguide potential, i.e., by creating a potential barrier 
of the form U{qi) = — \^y{qi) —^y{qi)^ • A natural case is that in 
which ^y{qi) either vanishes or is an irrelevant constant energy 
shift. U{qi) is then the potential required to flatten out the depth 
of the global potential of y. In addition, the acceleration of the guided 
matter-waves towards the region of high- curvature is prevented. To 
explore in detail this possibility, we consider an elliptical trap^°'^^'^^, 
associated with the curve r{u) = {a cos w, b sin w), with a > I? > 0, and 
circumference L. The CIP in an elliptical trap reads 



r(«)= 



8m {b^cos^u + a^ sin^u) 



(11) 



The eccentricity of an ellipse is defined by £ = [l — (^/a)^]^G[0,l] and 
can be used to quantify the deformation from a circle (for which a = 
b, 8 = 0). For a ring of radius a = b (y, with 8 = 0), the curvature is 
K{qi) = 1/a and the CIP becomes constant, and the ground state 
density profile is uniform along the arc length qi. For y with e > 0, the 
CIP comes into play and creates two attractive double wells, centered 
around the points with higher curvature qi = {0,1/2} {b < a) and 



with the minimum value 



8m ^4 



. The extent to which geometry- 
induced effects can be cancelled out by painting a barrier 
U{qi) = —^y{qi) is illustrated in Fig. 3. Such cancellation is 
effective as long as the consistency conditions for the dimensional 
reduction hold, which ceases to be the case as 8 is increased while the 
transverse width Gq remains fixed. 

The ground state density profile is a fairly robust quantity, but we 
note that this compensation is efficient as well for dynamical pro- 
cesses involving all spectral properties of the waveguide. Consider 



(a) 



(b) 



(c) 



Figure 4 | Curvature- induced suppression of temporal Talbot oscillations. Time evolution of the density profile n(qi,t) = j dqj^ n{qi,q^,t) of an initially 
tightly-localized wavepacket released in a two-dimensional elliptical waveguide, (a) For a ring trap (e = 0) /7(^i,t) exhibits Talbot oscillations as a result of 
the quadratic dispersion relation (left). Two Talbot oscillations are displayed, (b) Whenever e > 0, the CIP lifts the degeneracies in the spectrum and 
suppresses Talbot oscillations (e = 0.9). (c) The CIP can be cancelled out by modulating the depth of the trap (e = 0.9). L = ISOctq in all cases and 
hence the revival time Tr = mVl{%h) is constant for different values of e. 
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the time evolution of the density profile of an initially localized 
wavepacket released in the elliptical trap, displayed in Fig. 4. For a 
ring trap, where e = 0, the evolution of the density profile n{qi, t) 
weaves a highly structured interference pattern with "scars" in the 
plane (qi, t), known as a "quantum carpet Such quantum carpets 
exhibit a temporal analogue of the Talbot effect: in wave optics, the 
near-field diffraction pattern of a wave incident upon a periodic 
grating is characterized by a spatial periodicity^^'^^. The quantum 
dynamics of an initially localised wavepacket which is released in a 
two-dimensional ring trap exhibits a periodic revival of the initial 
state, with period t^^ = mVl{Tih) (the temporal analogue of optical 
Talbot oscillations). The reconstruction of the density profile at t = 0 
can be traced back to the quadratic dispersion relation of the trap 
and the degeneracies it entails^^'^°. This phenomenon has been 
experimentally observed in a variety of systems^^'^^. In a two-dimen- 
sional elliptical waveguide, the spectrum is modified and the disper- 
sion relation ceases to be quadratic. For g > 0, the geometry-induced 
potential lifts the degeneracy in the spectrum, leading to the sup- 
pression of Talbot oscillations. Nonetheless, the dynamics corres- 
ponding to a ring trap can be effectively recovered in an ellipcal 
trap with £ > 0 after compensating the depth of the waveguide 
potential. The reapperance of Talbot oscillations in compensated 
elliptical waveguides signals the isospectral properties with respect 
to the ring trap, illustrating the suppression of curvature- induced 
effects. 

Discussion 

The dynamics of matter waves in bent waveguides is severely dis- 
torted by the appearance of an attractive curvature-induced quantum 
potential. As matter wave circuits shrink in size and atomic velocities 
must be reduced to maintain single mode propagation, curvature- 
induced potentials impose practical limitations on minimum velo- 
cities, and methods to reduce their effects are needed. Using methods 
of supersymmetric quantum mechanics, we have introduced a 
framework to design sets of bent waveguides which share the same 
scattering properties. As a relevant example, an infinite family of 
reflectionless waveguides with a controllable number of bound states 
has been presented. As a complementary approach, we have dis- 
cussed the possibility of tailoring curvature-induced effects by con- 
trolling the depth of the waveguide trapping potential. Our 
discussion has been focused on the effects of curvature on guided 
matter-waves which are experimentally realizable by a variety of 
techniques including atom chip technology and the painted potential 
technique based on a time-averaged optical dipole potential^. Our 
results are however directly applicable to other systems such as 
optical waveguides and photonic lattices^^, in which curvature- 
induced potentials^^'^^, reflectionless potentials^^ and concepts of 
supersymmetric quantum mechanics^^'^^ have already been imple- 
mented in the laboratory, and that provide a natural alternative 
platform to experimentally explore the interplay between geometry 
and supersymmetry in quantum mechanics. 

Methods 

Let us consider the case in which the superpotential depends on a collective set of 
parameters Uq, O = O(^i;ao). The partner potentials ^ are shape- invariant if they 
are related by ^ ^ (^i; ^o) =^ _ ((lu (^\) + where the residual term R{ai) is 
independent of qi and Ui = f{ao) is a new set of parameters obatined form by the 
action of the function/(*)- By iteration, one can construct the series of Hamiltonians 
{Hk\k = 0, 1, ...} with Ho = H_ and Hi = H+, such that H, = Ho + Yl[=i ^(^^)' 
with flfc = f{ao), i.e., obtained by the/function iterated k times. It follows that the 
squared curvatures of the SUSY partner curves {ys- 1> Ys] related by a constant shift 
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